GROUP INVARIANT SOLUTIONS WITHOUT TRANSVERSALITY 
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Abstract. We present a generahzation of Lie's method for finding the group invariant solutions 
to a system of partial differential equations. Our generalization relaxes the standard transversality 
assumption and encompasses the common situation where the reduced differential equations for the 
group invariant solutions involve both fewer dependent and independent variables. The theoretical 
basis for our method is provided by a general existence theorem for the invariant sections, both 
local and global, of a bundle on which a finite dimensional Lie group acts. A simple and natural 
extension of our characterization of invariant sections leads to an intrinsic characterization of the 
reduced equations for the group invariant solutions for a system of differential equations. The char- 
acterization of both the invariant sections and the reduced equations are summarized schematically 
by the kinematic and dynamic reduction diagrams and are illustrated by a number of examples 
from fluid mechanics, harmonic maps, and general relativity. This work also provides the theoretical 
foundations for a further detailed study of the reduced equations for group invariant solutions. 

Keywords. Lie symmetry reduction, group invariant solutions, kinematic reduction diagram, dy- 
namic reduction diagram. 
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1. Introduction. Lie's metliod of symmetry reduction for finding the group invariant solutions to 
partial differential equations is widely recognized as one of the most general and effective methods for 
obtaining exact solutions of non-linear partial differential equations. In recent years Lie's method 
has been described in a number of excellent texts and survey articles (see, for example, Bluman 
and Kumci [10], Olvcr [29], Stcphani [36], Wintcrnitz [41]) and has been systematically applied to 
differential equations arising in a broad spectrum of disciplines (see, for example, Ibragimov [23] or 
Rogers and Shadwick [34]). It came, therefore, as quite a surprise to the present authors that Lie's 
method, as it is conventionally described, does not provide an appropriate theoretical framework 
for the derivation of such celebrated invariant solutions as the Schwarzschild solution of the vacuum 
Einstein equations, the instanton and monopole solutions in Yang-Mills theory or the Veronese map 
for the harmonic map equations. The primary objectives of this paper are to focus attention on 
this deficiency in the literature on Lie's method, to describe the elementary steps needed to correct 
this problem, and to give a precise formulation of the reduced differential equations for the group 
invariant solutions which arise from this generalization of Lie's method. 

A second impetus for the present article is to provide the foundations for a systematic study of 
the interplay between the formal geometric properties of a system of differential equations, such as 
the conservation laws, symmetries, Hamiltonian structures, variational principles, local solvability, 
formal integrability and so on, and those same properties of the reduced equations for the group 
invariant solutions. Two problems merit special attention. First, one can interpret the principle 
of symmetric criticality [32], [33] as the problem of determining those group actions for which the 
reduced equations of a system of Euler-Lagrange equations are derivable from a canonically defined 
Lagrangian. Our previous work [2] on this problem, and the closely related problem of reduction 
of conservation laws, was cast entirely within the context of transverse group actions. Therefore, 
in order to extend our results to include the reductions that one encounters in field theory and 
differential geometry, one needs the more general description of Lie symmetry reduction obtained 
here. Secondly, there do not appear to be any general theorems in the literature which insure the 
local existence of group invariant solutions to differential equations; however, as one step in this 
direction the results presented here can be used to determine when a system of differential equations 
of Cauchy-Kovalevskaya type remain of Cauchy-Kovalevskaya type under reduction [4]. 

We begin by quickly reviewing the salient steps of Lie's method and then comparing Lie's method 
with the standard derivation of the Schwarzschild solution of the vacuum Einstein equations. This 
will clearly demonstrate the difficulties with the classical Lie approach. In section 3 we describe, 
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in detail, a general method for characterizing the group invariant sections of a given bundle. In 
section 4 the reduced equations for the group invariant solutions are constructed in the case where 
reduction in both the number of independent and dependent variables can occur. We define the 
residual symmetry group of the reduced equations in section 5. In section 6 we illustrate, at some 
length, these results with a variety of examples. In the appendix we briefly outline some of the 
technical issues underlying the general theory of Lie symmetry reduction for the group invariant 
solutions of differential equations. 

2. Lie's Method for Group Invariant Solutions. Consider a system of second-order partial 
differential equations 

A^(x\k",<,<^-)-0 (2.1) 

for the m unknown functions m", a — 1, . . . , to, as functions of the n independent variables x% 
i ~ I, . . . , n. As usual, uf and uf, denote the first and second order partial derivatives of the 
functions u". We have assumed that the equations (2.1) are second-order and that the number 
of equations coincides with the number of unknown functions strictly for the sake of simplicity. A 
fundamental feature of Lie's entire approach to symmetry reduction of differential equations, and one 
that contributes greatly to its broad applicability, is that the Lie algebra of infinitesimal symmetries 
of a system of differential equations can be systematically and readily determined. We are not so 
much concerned with this aspect of Lie's work and accordingly assume that the symmetry algebra 
of (2.1) is given. Now let F be a finite dimensional Lie subalgebra of the symmetry algebra of (2.1), 
generated by vector fields 

^a-C(-^)|^+r;."(x^n^)|^, (2.2) 

where a = 1, . . . , p. A map s: R" -^ R™ given by u" = s"(x*) is said to be invariant under the 
Lie algebra F if the graph is invariant under the local flows of the vector fields (2.2). One finds this 
to be the case if and only if the functions s"(a;*) satisfy the infinitesimal invariance equations 

ci^n^^v:ix',s''{x^)) (2.3) 

for all a = 1, 2, . . . , p. The method of Lie symmetry reduction consists of explicitly solving the 
infinitesimal invariance equations (2.3) and substituting the solutions of (2.3) into (2.1) to derive 
the reduced equations for the F invariant solutions. 

In order to solve (2.3) it is customarily assumed (see, for example, Olver [29], Ovsiannikov [30], 
or Winternitz [41]) that the rank of the matrix [^^(2^"')] is constant, say q, and that the Lie algebra 
of vector fields satisfies the local transversality condition 
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k[C{x')]=rank[C{x'),rj:ix^,u'^)]. (2.4) 
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Granted (2.4), it then follows that there exist local coordinates 

F=F(x^'), x''^x^{x3) and w" = '^"(x^u^), (2.5) 

on the space of independent and dependent variables, where r == 1, . . . , n — g, fc = 1, . . . , q, and 
a = I, ... , TO, such that, in these new coordinates, the vector fields Va take the form 

K = E^l(^^-')|-7• (2-6) 

1=1 

The coordinate functions x^ and w" are the infinitesimal invariants for the Lie algebra of vector 
fields r. In these coordinates the infinitesimal invariance equations (2.3) for w" = w"(x'",x'^) can 
be explicitly integrated to give v" = w"(x''), where the w"(x'") are arbitrary smooth functions. One 
now inverts the relations (2.5) to find that the explicit solutions to (2.3) are given by 

s"(x'',i'=) =u"(x^x'=,w"(x'')). (2.7) 

Finally one substitutes (2.7) into the differential equations (2.1) to arrive at the reduced system of 
differential equations 

A/3(F,t;",<,O = 0. (2.8) 

Every solution of (2.8) therefore determines, by (2.7), a solution of (2.1) which also satisfies the 
invariance condition (2.3). In many applications of Lie reduction one picks the Lie algebra of vector 
fields (2.2) so that q = n — 1 in which case there is only one independent invariant i; on Af and (2.8) 
is a system of ordinary differential equations. 

For the vacuum Einstein equations the independent variables a;*, i ~ 0, . . . , 3, are the local 
coordinates on a 4 dimensional spacetime, the dependent variables are the 10 components gij of 
the spacetime metric and the differential equations (2.1) are given by the vanishing of the Einstein 

tensor G*-' = 0. In the case of the spherically symmetric, stationary solutions to the vacuum Einstein 

d 
equations the relevant infinitesimal symmetry generators on spacetime are Vq = -tt—^, 

ox" 

OT^ OX'^ OX^ OX"^ ox^ ox^ 

and the symmetry conditions, as represented by the Killing equations Cvadij = 0, lead to the familiar 
ansatz (in spherical coordinates) 

ds^ = A{r)df + B{r)dtdr + C{r)dr^ + D{r){d(p^ + s\n{(t>f dO^) . (2.9) 
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The substitution of (2.9) into the field equations leads to a system of ODE whose general solution 
leads to the Schwarzsehild solution to the vacuum Einstein field equations. 

What happens if we attempt to derive the Schwarzsehild solution using the classical Lie ansatz 
(2.7)? To begin, it is necessary to lift the vector fields Va to the space of independent and dependent 
variables in order to account for the induced action of the infinitesimal spacetime transformations 
on the components of the metric. These lifted vector fields are Vq — Vq and 

^ dV^ d 

14^14 -2^,,—. (2.10) 

In terms of these lifted vector fields, the infinitesimal invariance equations (2.3) then coincide exactly 
with the Killing equations. However, (2.7) cannot possibly coincide with (2.9) since the latter 
contains only 4 arbitrary functions A{r), B(r), C(r), D{r) whereas (2.7) would imply that the 
general stationary rotationally invariant metric depends upon 10 arbitrary functions of r. This 
discrepancy is easily accounted for — in this example 

rank [Vo,Vi,V2,V3] = 3 while rank [Vo , Vi , 1>2 , V3] - 4, 

and hence the local transversality condition (2.4) does not hold. Indeed, whenever the local transver- 
sality condition fails, the general solution to the infinitesimal invariance equation will depend upon 
fewer arbitrary functions than the original number of dependent variables. The reduced differential 
equations will be a system of equations with both fewer independent and dependent variables. 

We remark that in many of the exhaustive classifications of invariant solutions using Lie reduction 
either the number of independent variables is 2 and hence, typically, the number of vector fields Va 
is one, or there is just a single dependent variable and (2.1) is a scalar partial differential equation. 
In either circumstance the local transversality condition is normally satisfied and the ansatz (2.7) 
gives the correct solution to the infinitesimal invariance equation (2.3). However, once the number 
of independent and dependent variables exceed these minimal thresholds, as is the case in most 
physical field theories, the local transversality condition is likely to fail. 

3. An Existence Theorem for Invariant Sections. Let M be an n-dimensional manifold and 
TT : E -^ M a bundle over M . In our applications to Lie symmetry reduction the manifold M serves as 
the space of independent variables and the bundle E plays the role of the total space of independent 
and dependent variables. We refer to points of M with local coordinates (x*) and to points of E 
with local coordinates (a;*, u"), for which the projection map tt is given by 7r(a;*, u") — (a;*). In many 
applications E either is a trivial bundle E = M x N, a, vector bundle over M, or a fiber bundle over 
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M with finite dimensional structure group. However, for the purposes of this paper one need only 
suppose that tt is a smooth submersion. We let E^ = 7r^^(x) denote the fiber of E over the point 
X e M. 

Now let G be a finite dimensional Lie group which acts smoothly on E. We assume that G acts 
projectably on E in the sense that the action of each element of G is a fiber preserving transformation 
on E — if P, 9 lie in a common fiber, then so do g ■ p and g ■ q. Consequently, there is a smooth 
induced action of G on M. The action of G on the space of sections of E is then given by 

{g-s){x)^g-s{g-'-x). (3.1) 

for each smooth section s: M -^ E. 

A section s is invariant ii g ■ s ~ s for all g G G. More generally, we have the following definition. 

Definition 3.1. Let G be a smooth projectable group action on the bundle tt: E ^ M and let 
U d M be open. Then a smooth section s: U ^ E is G invariant, if for all x ^ U and g ^ G such 
that g ■ X ^ U, 

s{9-x) ^ g- s{x). (3.2) 

Let F be the Lie algebra of vector fields on E which are the infinitesimal generators for the action 
of G on E. Since the action of G is assumed projectable, any basis Va, a = 1, ... , p assumes the 
local coordinate form (2.2). If gt is a one-parameter subgroup of G with associated infinitesimal 
generator Va on E, then by differentiating the invariance condition s{gt ■ x) = gt ■ s{x) one finds that 
the component functions s"(x*) satisfy the infinitesimal invariance condition (2.3). If s is globally 
defined on all of M and if G is connected, then the infinitesimal invariance criterion (2.3) implies 
(3.2). This may not be true if G is not connected or if s is only defined on a proper open subset of 
M. 

For the purposes of finding group invariant solutions of differential equations, we shall take the 
group G to be a symmetry group of the given system of differential equations. The task at hand 
is to explicitly identify the space of G invariant sections of E with sections of an auxiliary bundle 
ti'kg • '*G (E) -^ M and to construct the differential equations for the G invariant sections as a reduced 
system of differential equations on the sections of tt^^ : k.q{E) — > M . 

Our characterization of the G invariant sections of E is based upon the following key observation. 
Suppose that p Cz E and that there is a G invariant section s : U ^ E with s{x) = p, where x e [/. 
Let Gx ^{gEG\g-x = x} he the isotropy subgroup of G at x. Then, for every g £ Gx, we 
compute 

g-p^ g- s{x) :^ s{g- x) ^ s{x) ^p. (3.3) 
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This equation shows that the isotropy subgroup Gx constrains the admissible values that an invariant 
section can assume at the point x. Accordingly, we define the kinematic bundle Kg{E) for the 
action of G on E by 



Ka{E) = y KaAE) 
xeM 

KgAE) = {p<^Ex\g-p = p for all geG^}. (3.4) 



xeM 
where 



It is easy to check that Kq (E) is a G invariant subset of E and therefore the action of G restricts to 
an action on Kg{E). 

Let M = M/G and Kg{E) — Ka{E)/G be the quotient spaces for the actions of G on M and 
Ko{E). We define the kinematic reduction diagram for the action of G on E to be the 
commutative diagram 

ka{E) ^^^^ Ka{E) — ^-^ E 

(3.5) 



1- 



M ^-^^^ M -^^ M. 
In this diagram t is the inclusion map of the kinematic bundle k.g{E) into i<^, id: M ^ M is the 
identity map, the maps q^ and q„<j are the projection maps to the quotient spaces and tt^^ is the 
surjective map induced by tt. The next lemma summarizes two of the key properties of the kinematic 
reduction diagram. 

Lemma 3.2. Let G act projectably on E. 

[i] Let p G Ka{E) and g G G. If Tr{g ■ p) = Tr{p), then g ■ p = p. 

[ii] If p (z Kg{E) and x G M satisfy Traip) — qM{x), then there is a unique point p G Kg{E) such 

that (^Koip) ^ P ^^d n{p) = X. 

Proof, [i] Let x = tt{p). If TT(g ■ p) = tt{p), then g ■ x = x and therefore, since p G Ko,^{E), we 
conclude that g ■ p ~ p. 

[ii] Since q^g : Kg{E) -^ Ka{E) is surjective, there is a point po G Hg{E) which projects to p. Let 
2^0 ~ '''(po)- Then (\m{xi:i) = (\m{x) and hence, by definition of the quotient map q^, there is a g G G 
such that g ■ xq — x. The point p ~ g ■ pa projects under q^^ to p and to x under tt so that the 
existence of the point p is established. Suppose pi and p2 are two points in Kq{E) which project 
to p and X under q^^ and tt respectively. Then pi and p2 belong to the same fiber Kg.^{E) and 
are related by a group element (/ G G, that is, g • Pi = P2- Since n{pi) = 7r(p2), it follows that 
^(5 ■ Pi) — ""(Pi)- Since pi G Kg,^{E), we infer from [i] that g ■ pi ~ Pi and therefore pi — p2- I 
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This simple lemma imm.cdiatcly im.plies that every local section s : U ^ Kg{E), where U is an 
open subset of M, uniquely determines a G-invariant section s: U ^ Ka(E), where U = q^/(C/), 
such that 

(\f^G{s(x)) ^ S{q^,r{x)). (3.6) 

To insure that this correspondence between the G invariant sections of E and the sections of kg (E) 
extends to a correspondence between smooth sections it suffices to insure that tt^c : Kg (E) -^ M is 
a smooth bundle. 

Theorem 3.3. (Existence Theorem for G invariant sections) Suppose that E admits a 
kinematic reduction diagram (3.5) such that Hg{E) is an imbedded subbundle of E, the quotient 
spaces M and ka{E) are smooth manifolds, and n^^ : Kg{E) -^ M is a bundle. 

Let U be any open set in M and let U ~ q^j (U). Then (3.6) defines a one-to-one correspondence 
between the G invariant smooth sections s: U ^ E and the smooth sections s: [/ — > Kg{E). 

We can describe the kinematic reduction diagram in local coordinates as follows. Since tt^q : 
Ka{E) ^ M is a bundle we begin with local coordinates tt^^ : (x^,u") -^ (x'^) for Ka{E), where 
r =: 1, . . . , diniM and a ranges from 1 to the fiber dimension of ka{E). Since qM : Af -^ M is a 
submersion, we can use the coordinates x*" as part of a local coordinate system (x*", x'') on M. Here 
fc = 1, . . . , dim A'/ — dim Af and, for fixed values of x*", the points (x'', x*^) all lie on a common G orbit. 
As a consequence of Lemma 3.2 [ii] one can prove that q^^ restricts to a diffcomorphism between 
the fibers of Kc{E) and Kg{E) and hence one can use (x'",x'^,w°) as a system of local coordinates 
on Kg{E). Finally, let {Sf\x^,u°') -^ (x'",x'^) be a system of local coordinates on E. Since Kg{E) is 
an imbedded sub-bundle of E, the inclusion map t: Hg{E) -^ E assumes the form 

t(x^x^i;") ^ (x^x^t"(x^x^^''^)), (3.7) 



where the rank of the Jacobian matrix 
reduction diagram (3.5) becomes 



dv'' 



is maximal. In these coordinates the kinematic G 



(x'',w'') < ""'' (x^x'=,w'') '-^ (x'',x'=,i"(x'',x'=,w'^)) 



(3.8) 



(x'-) ^^^^~ (x^x'=^) — ^^ (£^x'=). 

These coordinates are readily constructed in most applications. If u" ~ s'^{Sf') is a local section of 
kc{E), then the corresponding G invariant section of E is given by 



s"(r,x^) = i"(x^x^s"(x'')). (3.9) 
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Notice that when t is the identity map, (3.9) reduces to (2.7). The formula (3.9) is the full and 
proper generalization of the classical Lie prescription (2. 7) for infinitesimally invariant sections of 
transverse actions. 

In general the fiber dimension of Kg{E) will be less than that of E, while the fiber dimension of 
kc{E) is always the same as that of Kc{E). Thus, in our description of the G invariant sections of 
-E, fiber reduction, or reduction in the number of dependent variables, occurs in the right square 
of the diagram (3.5) while base reduction, or reduction in the number of independent variables, 
occurs in the left square of (3.5). 

We now consider the case of an infinitesimal group action on E, defined directly by a p-dimcnsional 
Lie algebra F of vector fields (2.2). These vector fields need not be the infinitesimal generators of a 
global action of a Lie group G on E. If the rank of the coefficient matrix [^^(a;-')] is q, then there 
are locally defined functions (/)°(x^), where e = 1, ... , p — q, such that 

^0:(x^)C(x-') = o. 

o=l 

Consequently, if we multiply the infinitesimal invariance equation (2.3) by the functions (f>'^{x^) and 
sum on a — 1, . . . ,p, we find that the invariant sections s"(x^) are constrained by the algebraic 
equations 

J2rA^')v:ix',s^{x')) = o. (3.10) 

a=l 

These conditions are the infinitesimal counterparts to equations (3.3) and accordingly we define the 
the infinitesimal kinematic bundle Kr{E) = [J^eM '^r.x{E), where 

KrAE) = { {x^uf") eE^Y. €{^')€{^\n^) - } 

a=l 

= {p<eE^\ Z{p) == for all Z e F such that 7r,(Z(p)) = }. (3.11) 

In most applications the algebraic conditions defining Hr{E) are easily solved. The Lie algebra 
of vector fields F restricts to a Lie algebra of vector fields on Kr {E) which now satisfies the infini- 
tesimal transversality condition (2.4). One then arrives at (3.8) as a local coordinate description of 
the infinitesimal kinematic diagram for F, where the coordinates (x^,w°) are now the infinitesimal 
invariants for the action of F on Kr(-E). 

It is not difficult to show that Kq^^E) C K,r.x{E), with equality holding whenever the isotropy 
group Gx is connected. 
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In the case where E is a vector bundle, the infinitesimal kinematic bundle appears in Fels and 
Olver [16]. For applications of the kinematic bundle to the classification of invariant tensors and 
spinors see [6] and [7]. 

4. Reduced Differential Equations for Group Invariant Solutions. Let G be a Lie group 
acting projectably on the bundle it: E ^ M and let A = be a system of G invariant differential 
equations for the sections of E. In order to describe geometrically the reduced equations A = 
for the G invariant solutions to A = we first formalize the definition of a system of differential 
equations. 

To this end, let tt'' : J^{E) ^ M be the fc-th order jet bundle oiir: E -^ M. A point a = f{s)(x) 
in J^{E) represents the values of a local section s and all its derivatives to order k at the point x £ M. 
Since G acts naturally on the space of sections of _E by (3.1), the action of G on i? can be lifted (or 
prolonged) to an action on J''(E) by setting 

9-^ ^ f{9- s){g- x), where a = f{s){x) 

Now let tt: I? -^ J^{E) be a vector bundle over J^{E) and suppose that the Lie group acts pro- 
jectably on I? in a manner which covers the action of G on J^{E). A differential operator is a 
section A : J^{E) -^ T). The differential operator A is G invariant if it is invariant in the sense of 
Definition 3.1, that is, 

g ■ A(a) = A(g • a) 

for all 5 e G and all points a E j'^ (E) . A section s oi E defined on an open set t/ C M is a solution 
to the differential equations A = if A{j''{s)(x)) = for all x £ U. 

Typically, the bundle T> -^ J^{E) is defined as the puUback bundle of a vector bundle V (on 
which G acts) over i? or M by the projections tt^ : J'^(E) -^ E or tt^ : J^{E) -^ M and the action 
of G on I? is the action jointly induced from J*'{E) and V . 

Our goal now is to construct a bundle T) -^ j''(ka(E)) and a differential operator A: J*'{ka{E)) 
-^ D such that the correspondence (3.6) restricts to a 1-1 correspondence between the G invariant 
solutions of A = and the solutions of A = 0. 

One might anticipate that the required bundle T) -^ J^{ko{E)) can be constructed by a direct 
application of kinematic reduction toT) ~* J^{E). However, one can readily check that the quotient 
space of J*^ {E) by the prolonged action of G does not in general coincide with the jet space J^ (ha (E) ) 
so that the kinematic reduction diagram for the action of G on 2? will not lead to a bundle over 
J''{kg{E)). For example, if G is the group acting on M x R ^- M by rotations in the base 
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M = R^ — {(0, 0)}, then J'^{E)/G is a 7 dimensional manifold whereas ,P{ka{E)) is 4 dimensional. 
This difficulty is easily circumvented by introducing the bundle of invariant k-jets 

lnv'^{E)={aeJHE)\a^f{s){xo), 

where s is a G invariant section defined in a neighborhood of xo }. 

(4.1) 

This bundle is studied in Olver [29] although the importance of these invariant jet spaces to the 
general theory of symmetry reduction of differential equations is not as widely acknowledged in the 
literature as it should be. 

The quotient space lnY''{E)/G coincides with the jet space J^{'ka{E)). We let I?i„^ —> lmr^{E) 
be the restriction of D to the bundle of invariant fc-jets and to this we now apply our reduction 
procedure to arrive at the dynamic reduction diagram 

^g(A„v) ^-^- Kc(A„0 -^-^ A„v -^^^^ V 

J^{ka{E)) ^^5i^ Inv'=(^) — ^^ lxw^{E) — ^ J''{E). 

Theorem 3.3 insures that there is a one-to-one correspondence between the G invariant sections of 
2?i„v -^ Inv {E) and the sections of KcC^^inv) -^ J^{ka{E)). Any G invariant differential operator 
A : J^ (E) -^ V restricts to a G invariant differential operator Ai„^ : Inv (E) -^ I?i„^ and thus 
determines a differential operator A: J^{ka{E)) -^ kaiT)). This is the reduced differential operator 
whose solutions describe the G invariant solutions for the original operator A. 

To describe diagram (4.2) in local coordinates, we begin with the coordinate description (3.8) of 
the kinematic reduction diagram and we let 

(z.r ~k a a a a a a \ 

denote the standard jet coordinates on J^{E). Since the invariant sections are parameterized by 
functions w" = v°-{3f), coordinates for Inv {E) are 

(^^x^«^ <,<„...) 

In accordance with (3.9), the inclusion map 

l: Itw^{E)^ J^{E) 
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is given by 

L[X ,X ,V ,V^,V^g, . . .) ^ (X ,X ,U ,Uj, ,Ui,,U^g,Uj,i^,Uf.i, . . .), (4-3) 

where by a formal application of the chain rule, 



u" 


c.,~r ^^ as a ^'^ ^'^ a a di" 


<s 


di'-dx' dv'^di'^ dv<^di^"' dv^dv''^'"' dv^ '''' 



and so on. The quotient map 

is given simply by 

Next let f '^ be a local frame field for the vector bundle V. The differential operator A : J''{E) -^ V 
can be written in terms of the standard coordinates on J^ (E) and in this local frame by 

A==A^(x^x^«^<,<,<„<fc,<^,...)f■". (4.4) 

The restriction of A to Inv (E) defines the section Ajnv : Inv (E) -^ I?i„^ by 

AI„v = AI„v,A(^^x^^.^^;^,l;^„...)f'', (4.5) 

where the component functions Ainv,A(a;'^,x'^,f",f",f"^, . . .) are defined as the composition of the 
maps (4.3) and the component maps A^. Since A is a G invariant differential operator, Ainv is a G 
invariant differential operator and hence necessarily factors through the kinematic bundle Kcl^^inv); 

Ai„v:Inv'=(£;)^«;G(A„v). 

Our general existence theory for invariant sections implies that we can also find a locally defined, 
G invariant frame f J^.^ for Kg(2?i„v)- The inclusion map Kq{'Di„^) -^ I?i„ is represented by writing 
each vector f JJ^ as a linear combination of the vectors f "^, 

where the coefficients Mj( are functions on Inv (E) . The invariant operator Ainv can be expressed 
as 

)(x^x^^;^^;,^<„...)ff„,. 
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Finally, the G invariant frame f J^^ determines a frame f^ on kc{E), the invariance of A implies 
that the component functions Ainv,Q are necessarily independent of the parametric variables x , 
that is, 

and the reduced differential operator is 

A = AQ(F,«^<,<,...)f^. 

At first sight, this general framework may appear to be rather cumbersome and overly compli- 
cated. However, as we shall see in examples, every square in the dynamic reduction diagram (4.2) 
actually corresponds to the individual steps that one performs in practice. 

5. The Automorphism Group of the Kinematic Bundle. Let (3 be the full group of pro- 
jectable symmetries on E for a given system of differential equations on J''{E) and let C? C 25 be 
a fixed subgroup for which the group invariant solutions are sought. It is commonly noted (again, 
within the context of reduction with transversality) that Nor(G, ©), the normalizer of G in ©, pre- 
serves the space of invariant sections and that Nor(G, >S)/G is a symmetry group of the reduced 
equations. However, because this is a purely algebraic construction which does not take into account 
the action of © on E, this construction may not yield the largest possible residual symmetry group 
or may result in a residual group which does not act effectively on ka{E). These difficulties are 
easily resolved. We let Op(G) denote orbit of G through a point p E E. 

Definition 5.1. Let & he a group of fiber-preserving transformations acting on tt: E ^ M and let 
G be a subgroup of &. Assume that E admits a kinematic reduction diagram (3.5) for the action of 
G onE. 

[i] The automorphism group © for the kinematic bundle ir : Kg (E) ^ M is the subgroup of (S 
which stabilizes the set of all the G orbits in Kg{E), that is, 

6 = {ae 0|a-Op(G) = Oa.p(G) and a'^ ' Op(G) = 0„-i.p(G) /or all p e Ka{E)}. 

(5.1) 

[ii] The global isotropy subgroup of(&, as it acts on the space of G orbits of Kg{E), is 

©* = { a e (5 I a • Op(G) = Op(G) for all p (^ Kg{E)] . (5.2) 

[iii] The residual symmetry group is ©cff ~ ©/©*. 

The key property of 0* is that it is the largest subgroup of © with exactly the same reduction 
diagram and invariant sections as G. This is an important interpretation of the group C5* — from 
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the viewpoint of kinematic reduction, one should generally replace the group G by the group 25*. 
For computational purposes, it is often advantageous to use the fact that ©* fixes every G invariant 
section of E. It is not difficult to check that Nor((5*, (S) = (S, that the quotient group ©off' = ©/©* 
acts effectively and projectably on the reduced bundle ka{E) -^ M and that, if (S is a symmetry 
group of of a differential operator A, then ©off is always a symmetry group of the reduced differential 
operator A. 

Similarly, if tj is a Lie algebra of projectable vector fields on E and F C ^, we define the 
infinitesimal automorphism algebra of Kr [E) as the Lie subalgebra of vector fields given by 

g = {Y (^g\[Z,Y]p(^iiY>(m{T){p) for all p e Kp(S) and all Z e F }, (5.3) 

and the associated isotropy subalgebra for Kp{E) 

g* ={Y e g \Yp e spa.n(T){p) ior all p e k^{E) } . (5.4) 

When © is a finite dimensional Lie group and g = F(©), then it is readily checked that g — F(©) 
and g* = F((S*). 

Since the automorphism group © acts on the fc-jets of invariant sections Inv q{E), this group also 
plays an important role in dynamic reduction. Specifically, let us suppose that (S acts on the vector 
bundle V -^ J^{E) and that A: J^{E) ^ V \s a & invariant section. Then Ai„,: Inv''(^) -^ Vj„^ 
is always invariant under the action of © and accordingly the operator Ai„^ always factors through 
the kinematic bundle for the action of © on I?i„v, where for a e Inv (E), 

«;e,,(A„v) = {Ae A„v,a|g-A = A for aU .9 € ©.}. 

We note that 

k8(A„v)c«;g (!?:„) 

and consequently one can refine the dynamic reduction diagram from (4.2) to 

4 4 4 b 

J''{ka{E)) ^-5i^^ lxw^{E) — ^^ lxw^{E) — ^ J^{E), 

where the quotient maps to the left are still by the action of G. 

Given the actions of G on tt : E ^ M and also © on I? ^ J*^ [E) , it sometimes happens that 

«0..(A„v) = 0. (5.5) 



Group Invariant Solutions without Transversality i4 

In this case every G invariant section of E is automatically a solution to A ~ for every © invariant 
operator A: J^{E) -^ V — such sections are called universal solutions. Previous work on this 
subject (see Bleecker [8], [9], Gaeta and Morando [19]) have emphasized a variational approach 
which, from the viewpoint of the dynamic reduction diagram and the automorphism group of the 
kinematic reduction diagram, may not always be necessary. 

6. Examples. In this section we find the kinematic and dynamic reduction diagrams for the group 
invariant solutions for some well-known systems of differential equations in applied mathematics, 
differential geometry, and mathematical physics. We begin by deriving the rotationally invariant 
solutions of the Eulcr equations for incompressible fluid flow. As noted by Olver [29] (p. 199), these 
solutions cannot be obtained by the classical Lie ansatz. The general theory of symmetry reduction 
without transversality leads to some interesting new classification problems for group invariant 
solutions which we briefly illustrate by presenting another reduction of the Euler equations. 

In our second set of examples we consider reductions of the harmonic map equations. We show 
the classic Veronese map from S"^ -^ S^ is an example of a universal solution. In Example 5.4 
we consider another symmetry reduction of the harmonic map equation which nicely illustrates the 
construction of the reduced kinematic space for quotient manifolds M with boundary. 

In our third set of examples, the Schwarzschild and plane wave solutions of the vacuum Einstein 
equations are re-examined in the context of symmetry reduction without transversality. We demon- 
strate the importance of the automorphism group in understanding the geometric properties of the 
kinematic bundle and, as well, qualitative features of the reduced equations. 

Finally, some elementary examples from mechanics are used to demonstrate the basic differences 
between symmetry reduction for group invariant solutions and symplectic reduction of Hamiltonian 
systems. 

Although space does not permit us to do so, the kinematic and dynamic reduction diagrams are 
also nicely illustrated by symmetry reduction of the Yang-Mills equations as found, for example, in 
[22], [25], [27]. In particular, it is interesting to note that the invariance properties of the classical 
instanton solution to the Yang-Mills equations (Jackiw and Rebbi [24]) imply that it is a universal 
solution in the sense of equation (5.5). 

Euler Equations for Incompressible Fluid Flow 

The Euler equations are a system of 4 first order equations in 4 independent and dependent 
variables. The underlying bundle E for these equations is the trivial bundle R^ x R** ^- R** with 
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coordinates (i,x, u,p) — > (i,x), where x — (x^ , x'^ , x^) and u = [u^ ,u'^,u^) and the equations are 

Ut + u-Vu = -Vp and V • u = 0. (6.1) 

The full symmetry group 25 of the Eulcr equations is well-known (see, for example, [23], [29], [34]) 

Example 6.1. Rotationally Invariant Solutions of the Euler Equations. The symme- 
try group of the Euler equations contains the group G ~ SO (3), acting on E, by 

R-{t,x,u,p) = {t,R-x,R-u,p) = {t,R]x\R]u^,p), (6.2) 

for R = (i?*) e SO(3). To insure that the action of G on the base R^ is regular we restrict to the 
open set M C R^ where ||x|| 7^ 0. The infinitesimal generators for this action are 

■d d 

Vk = £k^Jx'-Q- + £fe»j"'^- (6-3) 

We first construct the kinematic reduction diagram for this action. For a given point .tq = 
(to, xo) G M , the isotropy subgroup Gx^ for the action of G on M is the subgroup SO(2)xo C SO(3) 
which fixes the vector xq in R"'. Since the only vectors invariant under all rotations about a given 
axis of rotation are vectors along the axis of rotation, we deduce that for xo G M , 

KG..o(^) = {(io,xo,u,p)|i?-u = u forah i? G SO(2)xo } 
= { (io,xo, u,p) I u = Axo for some A G R}. 

The same conclusion can be obtained by infinitesimal considerations. Indeed, the infinitesimal 
isotropy vector field at xo for the action on M is 

Zj = X^SkijX -^-j 

and therefore, if (t,x, u,p) G k,o.x{E), we must have, by (3.11), 

This implies that x x u = and so u is parallel to x. 

Either way, we conclude that Kq{E) is a two dimensional trivial bundle {t, x. A, B) -^ (t, x), where 
the inclusion map l : Kq (E) -^ E is 

L{t,x,A,B) — {t,x,u,p), where u = Ax and p = B. 
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The invariants for the action of G on M are t and r = \/ x^ + y^ + z^ so that the kinematic reduction 
diagram for the action of SO(3) on E is 



{t,r,A,B) < 


I^G 


- (i,x,Ai3) — 


— > {t,x,u,p) 


1-Sg 




TT 


7T 


(i,r) < 


(\M 


- (i,x) — 


-^ (i,x). 



(6.4) 



In accordance with equation (3.9), each section A = A{t,r) and B = B{t,r) of Ka{E) determines 
the rotationaUy invariant section 



u = A{r,t)x and p — B{r,t) 



(6.5) 



oi E. 



The computation of the reduced equations for the rotationaUy invariant solutions to the Euler 
equations now proceeds as foUows. From (6.5) we compute 



ui ^ Atx\ u) = AS] + A, 
so that the Euler equations (6.1) become 

Atx' + Ax' {A5] + Ar^^) ^ - 
which simplify to the differential equations 

At + A{A + rAr) 



Br 



and Pi — Br — 



-Br— and 3A + rAr = 
r 



and 'iA + rAr = 



(6.6) 



(6.7) 



(6.8) 



on J {kQ{E)). These equations are readily integrated to give 

a a a? 
A^ ^ and B ^ 7r^+b 

for arbitrary functions a(t) and b{t) and the rotationaUy invariant solutions to the Euler equations 
are 



(6.9) 



a a a 
u = — :7X and p = j 



b. 



We note that for the Lie algebra of vector fields (6.3), the matrix on the right side of (2.4), namely 



r 


~x^ 


x" 





-u^ 


u^ 1 


a;3 





-x^ 


u^ 





-u^ 


-x^ 


x^ 





~m2 


m1 
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has full rank 3 whereas the matrix on the left side of (2.4), consisting of the first three columns of 
(2.4), has rank 2. The local transversality condition (2.4) fails and the solution (6.9) to the Eulcr 
equations cannot be obtained using the classical Lie prescription. 

To describe the derivation of the reduced equations in the context of invariant differential operators 

and the dynamic reduction diagram we introduce the bundle V = J^(E) x R^ x R with sections 

d 

— — (8) dt and dt and define the differential operator A on I? by 

A = [mJ + u'^ul + S'^pj] |— ® di + [ui\ dt. (6.10) 

This operator is invariant under the full symmetry group of the Euler equations. The induced action 
of G = SO(3) on J^E) is given by 

R ■ {t, x\u\p, u),pj) = (i, i?;x^ i?>^p, i?j,i?j<, i?>,.) where R e SO(3). 

Coordinates for the bundle of invariant jets Inv^(ii^) are (i, a;% A, At, Ar, B, Bt, Br) and (6.6) defines 
the inclusion map i : Inv [E) -^ J^{E). A basis for the G invariant sections of Vinv -^ Inv (E) is 
given by 

{^=x'^(g,dt and i^^dt. 

Let f ^ and f ^ be the corresponding sections of Kg(2?i„v)- 

We are now ready to work though the dynamic reduction diagram (4.2), starting with the Euler 
operator as a section A: J^{E) -^ V. Restricted to the invariant jet bundle Inv^(_E), A becomes 

A,„, = [Atx' + Ax^{A5] + Ar^) + B,^] |- ® di + [Sl{A6) + Ar^)] dt. 

Restricting A to Inv (E) is precisely the first step one takes in practice in computing the reduced 
equations and corresponds to the right most square in the dynamic reduction diagram. 

Next, because Ai„^ is G invariant it is necessarily a linear combination of the two invariant sections 
f ^ and f^ and therefore factors though the kinematic bundle Kcil^inv)- This means we can write 
Ainv as a section of Kg(2?i„v), namely, 

A:„. = [At + A{A + Ar^)+^Br] f ^ + [3^1 + AriSi^)] f 2 

This corresponds to the center commutative square in the dynamic reduction diagram (4.2) and 
coincides with the fact that equation (6.7) contained a common factor x* - a common factor which 
insures that the time evolution equation for u reduces to a single time evolution equation for A. 
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Finally, as a G invariant section of Kcl^^inv), a- bundle on which G always acts transversally, 
we are assured that Ainv descends to a differential operator A on the bundle J^{ko{E)). In this 
example this implies that the independent variables (t, x^) appear only though the invariants for the 
action of G on M, in this case t and r, and so 



A = [^t + A{A + rAr 



^]i' + [3A + rAr]~f^. 
r 



Example 6.2. A new reduction of the Euler equations. It is possible to give a complete 
classification of all possible symmetry reductions of the Euler equations (6.1) to a system of ordinary 
differential equations in three or fewer dependent variables [17]. A number of authors have obtained 
complete lists of reductions of various differential equations (see, for example, [i4],[i8], [21], [42]) but 
this particular classification of reductions of the Euler equations may be the first such classification 
of group invariant solutions which explicitly requires non-trivial isotropy in the group action on the 
space of independent variables. There are too many cases to list the results of this classification 
here, but we do present one more reduction of the Euler equations, one which does not seem to 
appear elsewhere in the literature. 

For this example it will be convenient to write x = {x,y,z) and u ~ {u,v,w). The infinitesimal 
generators for the group action are F = { Vq, Vi, V2 = Vx.a + Vy.p, V3 — Vy^a — Vx.p }, where 



Vo — xdx + ydy + zdz + udu + vdy + wdw + 2p9p, Vi = ydx — xdy + vdu — ud^, 

VxM = adx + adu - xadp, and Vyjj = (3dy + pdy - yPdp, 

Here a = a{t) and /3 — j3{t) are such that a/3 — d/3 — 0, or equivalently, 

a/3 — j3a — c ~ constant. (6-11) 

This condition insures that [V2 , V3] = so that F is indeed a finite dimensional Fie algebra of vector 
fields. In order that F have constant rank on the base space, we assume that xya 7^ or yz(3 ^ 0. 
The horizontal components of V2 and V3 are given by 





= 


■ a 13' 
-(i a 




'dx' 

dy_ 



SO that 



dx 
dy 



a -P 
(3 a 






where (5 = a^ + /3^, and therefore at the point (io, xq), the horizontal components of the vector field 



^ ,. a{to)V2-l3{h)Vz , l3{to)V2+a{to)V^ 
Z = Vi- yo ^j—^ h xq- 



6{to) 



Sito) 



(6.12) 
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vanish. The isotropy condition (3.11) defining the fiber of the kinematic bundle Kr^(i?) leads, from 
the coefficients of 9„, d^ and dp, to the relations 

ya~xl3 xa + y(3 ■ xa + yf3 x(3 - ya ■ 

v = r a^ p and u = a -\ fj. (6.13) 

00 00 

We therefore conclude that the kinematic bundle has fiber dimension 2 with fiber coordinates w and 
p. However, these coordinates are not invariant under the action of F on k^{E) and cannot be used 
in the local coordinate description (3.8) of the kinematic reduction diagram. Restricted to Ki-{E), 
the vector fields Vi become 

Vq = xdx + ydy + zd^ + wdw + 2pdp, V( = yd^ - xdy, 

V2 = adx + (idy — {ax + (3y)dp, and V^ = —(3dx + ady — (—/3a; + ay)dp. 

Note that these restricted vector fields now satisfy the infinitesimal transversality condition (2.4). 
Invariants for this action are t, 

A=- and B={2p+'^^^-^{x^ + y^))/z\ (6.14) 

z 

To verify that B satisfies V^iB) — V^{B) — one must use a/3 = a/3. The kinematic reduction 
diagram for the action of F on i*^ is therefore 

{t,A,B) < """^ {t,x,A,B) — ' — > (i,x,u,p) 

4 4 h 

(i) < (i,x) > (i,x), 

(\M id 

where the inclusion map t is defined by (6.13) and the solutions to (6.14) for w and p. The general 
invariant section is then, on putting a = ln6, 



u^x--y-, v^x-+y-, 

W = zA{tl p ^ _^^±^ (^2 ^ y2) ^ ^5(,)^ 

Note that the u and v components are uniquely determined from the isotropy conditions (6.12) and 
(6.13) and that the arbitrary functions Ait) and B{t) defining these invariant sections appear only 
in the w and p components. 
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We now turn to the dynamic reduction diagram. Since we are treating the Euler equations as 
the section (6.10) of the tensor bundle T) we can anticipate the form of Ai„ by computing the F 
invariant tensors of the form 

T ^ P du (E) dt + Q dy (g) dt + Rd^ (g) dt + Sdt. (6.15) 

The isotropy condition CzT = at xq, where Z is defined by (6.12), shows immediately that 
P = Q = from which it follows that 

f^ = z-— (g)di and i^ ^ dt 
ow 

are a basis for the F invariant fields of the type (6.15). This calculation shows that the 9„ (E) dt and 
dy (g) dt components of the reduced Euler equations must vanish identically and, consistent with this 
conclusion, one readily computes 

r(T ,&.^ ,C-.2 (ofci + /9/3) n . d , d , 1 / ; ,9 \„i / ,\„9 

A:„ - [2 + ( 2 )' - {-^) / ] [^^ ®dt + y-® dt] + {A + A^ + i3)f 1 + {a + A)f ^ 

= (i + A^ + B)f 1 + (ct + yl)f2. 

Thus, the reduced differential equations are 

A + A^ + B = and a + A = {) 

which determine A and B algebraically. In conclusion, for each choice of a and [3 there is precisely 
one F invariant solution to the Euler equations given by 

aa + PP a[3 — d/3 a/3 — 6l[3 aa + /3/3 aa + P(3 

" = ^ ^2 + ^2 - y ^2 + ^2 ' " = '"^ ^2 + /32 + y C.2 + /32 ' "" = "^^ ^2 + ^2 ' 

1 2 2\Q^'3! + /3/3 2('^^^PP ,aj3^l3a-.2 /aa + j3f3-.2 

P=-^[^ +y ) ^2^P2 + ^ 1^ a2 + ^2 + I ^2 + ^2 J "•^ia2+/32^ 

Harmonic Maps 

For our next examples we look at two well-known reductions of the harmonic map equation for 
maps between spheres. For these examples the bundle _E is S*" x S*™ -^ 5" which we realize as a 
subset of R"+i X R™+i by 

^ = { (x, u) e R"+^ X R^+i |x-x = u-u = l}. 
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Let G be a Lie subgroup of SO(n + 1), let p: G ^ SO(?Ti + 1) be a Lie group homomorphisni and 
define the action of G on _E by 

R ■ (x, u)^{R- X, p{R) ■ u) for R e G. 

The kinematic bundle for the G invariant sections of E has fiber 

Kg,„{E) = {(x, u) e i? I p{R) • u = u for aU i? e G such that i? • x = x }. 

We identify the jet space J'^(E) with a submanifold of j2(R"+i, R^+i) by 

J^(E) = { (x, u, ^^u, diju) e j2(R"+i, R"+i) | x-x = 1, u-u = 1, u-a,u = 0, u■^iJU+^^u■^JU = }. 

Since the harmonic map operator (or tension field) is a tangent vector to the target sphere 5™ at 
each point a e J^{E), we let 

V = {(a, A) e J^{E) X R"+i I u • A = 0}. (6.16) 

By combining Proposition LI. 17 (p. 19) and Lemma VIL1.2 (p. 129) in Eells and Ratto [15], it follows 
that one can write the harmonic map operator A : J^ (E) ^ I? as the map 

A(cr) = [A^"^'u" + x'x'u^j + nx'uf - Au"] — ^, (6.17) 

where 

A = Sc^plS'^ufu'^ - x'x^ufu'^] and A^"^'u" = -S'^uf^. 

This operator is invariant under the induced action of © = SO(n + 1) x SO(?7i + 1) on E. 

Example 6.3. Harmonic maps from 5*^ to 5"'. For our first example we take E = S'^ x S'^ ^ S"^ 
and we look for harmonic maps which are invariant under the standard action of SO (3) acting 
on S^. It can be proved that, up to conjugation, there are three distinct group homomorphisms 
p: SO(3) -^ SO(5), which lead to the following three possibilities for the infinitesimal generators of 
SO (5) acting on E. 

Vi = zdy -yd,, . Vi^ zdy - yd, - u^d^s + u^^d^^ , 

Case I \y2=xd,- zd^, Case II < ^2 = xd, - zd^ - u^d^i + u^d^^ , 

V^ = yd^ -xdy. [vs^ yd^ - xdy - u^d^^ + u^d^i ■ 

Vi = zdy — ydz + u^d^i — u^du2 + {u^ — \/3u^)du3 — u'^du* + VSu^duH, 
Case III { V2 = xd, - zd^ - u^d^i + {u'^ + V3u^)du2 + u^d^^ - u^d^4 - VSu^d^5 , 
V3 = ydx — xdy — 2u^du^ + u'^9„2 — u^9„3 + 2u^9„4. 
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In Case I, the map p is the constant map, and in Case II, p is the standard inchision of SO(3) into 
SO (5). The origin of the map p in Case III will be discussed shortly. 

Since SO(3) acts transitively on 5^, the orbit manifold M consists of a single point, the space 
of invariant sections is a finite dimensional manifold, and the reduced differential equations are 
algebraic equations. The kinematic bundles Hg{E) are determined in each case from the isotropy 
constraint 

xVi + yV2 + zV^ = 0. 

In Case I the action is transverse, the isotropy constraint is vacuous and the kinematic bundle is 
kg{E) = S'^ X S*^. The invariant sections are given by 

^,ix,y,z)^{A,B,C,D,E), 

where A, . . . , E are constants and A^ + B^ + C^ + D^ + E^ = 1. In Case II the kinematic bundle 
is 5*^ X S'^ and the invariant sections are 

$„(x, y, z) = {Ax, Ay, Az, B, C), 

where A, B, C are constants such that A^ + B^ + C^ = 1. We take A ^ 0, since otherwise $77 
becomes a special case of $/. In Case III, Ka(E) = 5*^ x { ±1 } and the invariant sections are 

1 /3 

<^iji{x,y,z) = AV3{xy,xz,yz,-{x^ -y^),— (x^ + y^ -22^)), 

where A = ±1. 

Direct substitution into (6.17) easily shows that the maps $j and ^j^ automatically satisfy the 
harmonic map equation. The map $jj is harmonic if and only if B = C = in which case $jj 
is either the identity map or the antipodal map on S*^ followed by the standard inclusion into S^. 
Despite the simplicity of these conclusions, it is nevertheless instructive to look at the corresponding 
dynamic reduction diagrams. 

In Case I, the invariant sections are constant and so 

liw^{E) = {(x, A) e R-'' X R^ I X • X = A • A = 1, } 
and 

Vinv = {{a, A) e Inv^{E) x R^ | A • A = 0}. 

The automorphism group for the kinematic bundle in this case is (S = SO (3) x SO (5) which acts 
on Vinv by 

{R, S) ■ (x, A, A) = (i? • X, 5 • A, S* • A) for R e SO(3) and S e SO(5). 
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The isotropy constraint for k^ (I'inv) forces A to be a multiple of A. Hence, by the tangency 
condition A • A = 0, we have A = and n^ ^l-^inv) — 0- This shows that the map (f>j is harmonic 
by symmetry considerations alone and moreover that it is a universal solution for any operator 
A: J'=(52 X S^)->V with SO(3) x SO(5) symmetry. 

In Case II, the harmonic map equations force _B = C = so that the maps ^n are not universal. 
Interestingly however, the standard and antipodal inclusions 5*^ -^ S* have a larger symmetry 
group, namely SO(3) x SO(2) C 25 and it is easily seen, using these larger symmetry groups, that 
the standard and antipodal inclusions are universal. It is a common phenomenon that the group 
invariant solutions to a system of differential equations possess a larger symmetry group than the 
original group used in their construction. 

In Case III one finds immediately that k^ „(I?inv) ~ and ^m is universal, again for any operator 
A: 7*^(52 X 5"') ^ P with SO(3) x SO(5) symmetry. 

The map ^ju is the classic Veronese map. The symmetry group defining it is based on a stan- 
dard irreducible representation of SO(3) which readily generalizes to give harmonic maps between 
various spheres of higher dimension. Specifically, starting with the standard action of SO(n) on 
V = R", consider the induced action on Sym^^(y), the space of rank k symmetric, trace- free ten- 
sors or, equivalently, on the space W ~ Ti.'^iV) of harmonic polynomials of degree k on V. The 
standard metric on W is invariant under this action of SO(n) and in this way one obtains a Lie 
group monomorphism p: SO(n) -^ SO(iV), where A^ = dim(Vl^) = (" ^.^ ) ^ 1- -'^'-"' t^xample, the 
polynomials 

u^^xy, ,u^^xz, u^^yz, u"^ = l/2{x^ - y^), u^ = Vi/6{x^ + y'^ - 2z^) 

form an orthogonal basis for 7i^(R^) and the action of SO(3) on this space determines the action 
of SO(3) on R"^ X R''' in Case III. For further examples see EeUs and Ratto [15] and Toth [38]. | 

Example 6.4. Harmonic Maps from 5*" to 5*". A basic result of Smith [35] states that 
each element of 7r„(S'") = Z can be represented by a harmonic map (with respect to the standard 
metric) provided n < 7 or n = 9. This result, which can be established by symmetry reduction of 
the harmonic map equation (see Eells and Ratto [15] and Urakawa [39]), illustrates a number of 
interesting features. First, we see that much of the general theory which wc have outlined could be 
extended to the case where M is a manifold with boundary and where the fibers of Ka{E) change 
topological type on the boundary. Secondly, we find that the invariant sections for the standard 
action of G = SO(n— 1) x SO(2) C SO(n-|- 1) on S'" are slightly more general than those considered 
in [15] and [39]. However, a simple analysis of the reduced equations, based upon Noether's theorem, 
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shows that the only sohitions to the reduced equations are essentiaUy those provided by the ansatz 
used by Eells and Ratto and Urakawa. 

If {R, S)eG = SO(n - 1) X SO(2) C SO(n + 1) and 

(x, y,u,v)eEc (R""^ X R2) X (R"-^ x R^), 

where ||x|p + ||y^|| = 1 and ||u|p + ||v|p = 1, then the action of G on _E = 5" x S*" is given by 



(i?,5)(x,y,u,v)-( 



The invariants for the action of G on the base R"+^ are r = ||x|| and s = ||y|| which, for points 
(x,y) e S", are related by r^ + s^ = 1, where r > and s > 0. The quotient manifold M — S'^/G 
is therefore diffeomorphic to the closed interval [0, 7r/2]. 

To describe the kinematic bundle Kq (E) we must consider separately those points in M for which 
(i) s = 0, (ii) s 7^ and r 7^ and (iii) r = 0, corresponding the left-hand boundary point, the 
interior points and the right-hand boundary points of M. For (x,0) G S*", the isotropy subalgebra 
is SO(n — l)x X SO(2) and the fiber of the kinematic bundle consists of a pair of points 

Kg,(>c,o)(£^) = {(x, 0, u,v)|u = ±x, and v = 0}. 

For points (x,y) G S" with r 7^ and s 7^ the isotropy group is SO(n — l)x x {/} and the fiber 
of the kinematic bundle is the ellipsoid of revolution 

KG,(.,y)(-B) = {(x,y,u,v)|u = Ax, where r^A^ + ||v||2 = 1}. 

X • u y ■ V y • V 

Invariant coordinates on Ko,(x,y){E) are A ~ — k—, B ~ — ^^ ^^d C* = 2 — : where y = 



(0,-y^y^), subject to 

The inclusion map from KG,(jc,y)(-E') to -E(x,y) is 

u = Ax and v = By + Cy^ . 



(6.18) 



At the points (0, y), the isotropy subalgebra is SO(n) x { / } and the fiber of the kinematic bundle 
is the circle 

«G.(o,.,(^)-{(0,y,u,v)|u = and ||v|| = l}. 
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The G invariant sections are therefore described, as maps $: R"+^ -^ R"+^, by 

$(x, y) = A(t)x + B(t)y + C(t)y^, (6.19) 

f s 
where t is the smooth function of (x, y) defined by cos(t) = -^ ^ and sin(i) = -^ -, and where 

cos'^ {t)A'^{t) + sin^ {t){B^{t) +C^{t)) = 1. The isotropy conditions at the boundary of M imply that 
the functions A, B and C are subject to the boundary conditions 

A{0)=±1, B{0)=0, C(0)-0, and A(|) - 0, B{^r + 0(^^ = 1. (6.20) 

The invariant sections considered in [15] and [39] correspond to C(t) = 0. Note that the space of 
invariant sections (6.19) is preserved by rotations in the v plane, that is, rotations in the BC plane 
and therefore SO (2) C ©cff. 

By computing kg(I?i„^) we deduce that the restricted harmonic operator A,„^ is of the form 

A,„. = A. (x- 1-) + A, (y i-) + Ac (y^ • 1^), 
where the tangency condition (6.16) reduces to 

r^AA^ + s^BAs + s^CAc = 0. 

A series of straightforward calculations, using (6.17), now shows that the coefficients of the reduced 

operator A are 

- •• / sin(t) cos(i)s ; . < , 

^A-^-A+ (n )-( - ^^)A + nA- \A, 

cos(i) sm(t) ' 

~ - , , ^ sin(t) cos(t) s • , ^ 

A^^-B+{{n-2)^^-3^^)B + nB-XB, 6.21 

cos(t) sm(t) 

Ae = -C + (in - 2)^ - 3^)C + nC - AC, 
^ cos(t) sm(t) ' 

where 

A = cos^ (t)i2 + sin^ (t) {B^ + C^)+2 cos(i) sin(t) {-AA + BB + CC) 

+ (n - 1)^2 + 2(^2 + c^) _ cos2(i)A2 - sin2(t)(B2 + C^). 

To analyze these equations, we first invoke the principle of symmetric criticality and the formulas 
in [2] for the reduced Lagrangian to conclude that these equations are the Euler-Lagrange equations 
for the reduced Lagrangian 

L= -cos(t)"^2sin(t)A(it 
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subject, of course, to the constraint (6.18). From knowledge of the automorphism group of the 
kinematic bundle we know that this Lagrangian is invariant under rotations in the BC plane and 
this leads to the first integral 

J = cos(t)"-2 sin(i)3(BC' - CB) 

for (6.21). By the boundary conditions (6.20), J must vanish identically. Thus C{t) ~ fJ-B{t), for 
some constant /i and therefore a rotation in the v, v^'^ plane will rotate the general invariant section 
(6.19) into the section with C{t) = 0. We then have r'^A'^ + s^B^ = 1 and the change of variables 

^ ' cos(i) ^ ' sin(f) 

converts the reduced operator (6.21) into the form found in [15] or [39]. | 

General Relativity 

We now turn to some examples of Lie symmetry reduction in general relativity which we again 
examine from the viewpoint of the kinematic and dynamic reduction diagrams. To study reductions 
of the Einstein field equations, we take the bundle E to be the bundle Q{M) of quadratic forms, 
with Lorentz signature, on a 4-dimensional manifold M . A section of E then corresponds to a choice 
of Lorentz metric on M. We view the Einstein tensor 

d d 

A = G*-' {ghk , ghk,i , 9hkdm)-g-i «) -^ 

formally as a section oiV ^ J'^(E), where V is puUback of y == Sym^(T(M)) to the bundle of 2-jets 
J'^{E). The operator A is invariant under the Lie pseudo-group (S of all local diffeomorphisms of 
M. 

Let Divg be the covariant divergence operator (defined by the metric connection for g) acting on 
(1,1) tensors, 

Divg(S) = WiS}dx\ 

The contracted Bianchi identity is Div^ A^ = 0, where A'' is the operator obtained from A by 
lowering an index with the metric. 

The first point we wish to underscore with the following examples is that the kinematic reduction 
diagram gives a remarkably efficient means of solving the Killing equations for the determination 
of the invariant metrics. Secondly, we show that discrete symmetries, which will not change the 
dimension of the reduced spacetime M, can lead to isotropy constraints which reduce the fiber 
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dimension of the kinematic bundle. Thirdly, for G invariant metrics, the divergence operator Divg 
is a G invariant operator to which the dynamical reduction procedure can be applied to obtain the 
reduction of the contracted Bianchi identities for the reduced equations. Throughout, we emphasize 
the importance of the residual symmetry group in analyzing the reductions of the field equations. 

Finally, we remark that our conclusions in these examples are not restricted to the Einstein 
equations but in fact hold for any generally covariant metric field theories derivable from a variational 
principle. 

Example 6.5. Spherically Symmetric and Stationary, Spherically Symmetric Reduc- 
tions . We begin by looking at spherically symmetric solutions on the four dimensional manifold 
M = R X (R^ — { }), with coordinates (x*) = (t, x, y, z) for i = 0, 1, 2, 3. Although this is a very 
well-understood example, it is nevertheless instructive to consider it within the general theory of Lie 
symmetry reduction of differential equations. The infinitesimal generators for G = SO(3) are given 
by (2.10) and, just as in Example 6.1, we find that the infinitesimal isotropy constraint defining 

Kg,x(-E) = Kr,x(-E) is 



sokijX gii 



d 
dgij 



0, 



or, in terms of matrices, 



where 



ga + a*g = 0, 



(6.22) 
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lations are easily solved to give 
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(6.23) 



The fiber of the kinematic bundle Ko.x{E) is therefore parameterized by four variables A, B, C, D. 
Since these variables are invariants for the action of G restricted to Kq (E) and since the invariants 
for the action of SO(3) on M are t and r, the kinematic reduction diagram for the action of SO(3) 
on the bundle of Lorentz metrics is 

(In 



(t,r,A,B,C,D) 



{t,r) 



(IM 



{x\A,B,C,D) 



(x^) 



ix\g^J) 



(6.24) 



id 



(^0, 
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where the inclusion map l is given by (6.23). 

Consequently, the most general rotationally invariant metric on M is 

ds'^ = A{t,r)dt'^ + 2B{t,r)dt{x dx + y dy + z dz) + C{t,r){x dx + y dy + z dz)'^ + D{t,r){dx^ +dy'^ + dz'^). 

In standard spherical coordinates x ^ r cos 9 sin (/), y = r sin 9 sin (p, z = cos (p this takes the familiar 
form (on re-defining the coefficients _B, C and D) 

ds^ = A{t, r)df + B{t, r)dtdr + C(i, r)dr^ + D{t, r)dVl? (6.25) 

where dVi^ ^ dcf)"^ + s\n^ (j)d9^. 

d 
If we enlarge the symmetry group to include time translations Vq = -—^ then the kinematic 

at 
reduction diagram becomes 

{r,A,B,C,D) ^^^^ {x\A,B,C,D) -^-^ {x\g,,) 



(6.26) 



(r) ^^^^ (x*) -i^ (x^. 



At first glance there appears to be little difference between the two diagrams (6.24) and (6.26), 
but a computation of the automorphism groups reveals a dramatic difference in the geometry of the 
reduced bundles ka{E) in (6.24) and (6.26). This difference is best explained in terms of general 
results on Kaluza-Klein reductions of metric theories as in, for example, Coquereaux and Jadczyk 
[13]. From our perspective, these authors show that when the action of G on M is simple in the 
sense that the isotropy groups Gx can all be conjugated in G to a fixed isotropy group Gx^ , then 
the reduced bundle k.q{E) is a product of three bundles over M, 

ka{E)^Q{M)®A{M)®Qin4K). (6.27) 

Here 

[i] Q{M) is the bundle of metrics on M. 

[ii] A{M) — A^(M) (g) (P X/f f)) , where P is the principal H bundle defined as the set of points in 

M with isotropy group Gxo'And H — Nor(Gj;„, G)/Gxo- 

[iii] Qiyiv{K) is the trivial bundle whose fiber consist of the G invariant metrics on the homogeneous 

space K ~ GJGxq- 

For (6.24) one computes the residual symmetry group ©off to be the diffeomorphism group of 
M = R X R+ and one finds that the coefficients A, i?, G transform as the components of a metric 
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on M and that D is a scalar field (which one identifies as a map into the space of SO (3) invariant 
metrics on S*^). Thus, for (6.24), wc find that 

Ka{E) = Q{M)®n, 

where 3? is a trivial line bundle over M. By contrast, for the diagram (6.26) the automorphism 
group © acts on M by 

r^f{r) and t^et + g{r), (6.28) 

where / G Diff(R+), g G C°°(R) and e e R*. Without going further into the details of the 
decomposition (6.27), we simply note that the variable t is now the fiber coordinate on the principle 
bundle P and that under the transformations (6.28) the coefficients of the metric (6.25), which are 
now functions of r alone, transform according to 

A{r) ^ e^A{f{r)), B{r) ^ 6[/'B(/(r)) + 2g'A{f{r))] 

C{r) ^ {ffC{f{r)) + f'g'B{f{r)) + {g'fA{I{r)) D{r) ^ D{f{r)). 

Consequently, the sections of ka (E) can be written as 

s{r) = [g{r),uj{r),h{r)], 

where 

g{r) = [C{r) - ^^] dr^, u{r) = ^j^dr (g, ?-, and h{r) ^ A{r)dt^ + D{r) dVL^ . 

Here g{r) is a metric on M, a;(r) is a connection on P pulled back to M, and h(r) is a map from M 
into the G invariant metrics on R x S'^. 

The detailed expression for the reduced operator A for the stationary, rotationally invariant 
metrics can be found in any introductory text on general relativity. Here we simply point out that 
by computing the action of G on Sym^ (TM), we can deduce that the reduced operator will have 
the form 

X xtfd d ~.rf,d d d d . xrr-t^ 9 ~.o,d d Id d . 

A = A — (g) h A ( — (g) 1 (g) — ) + A — ® h A ( — ® 1 ^ (g) — ), 

dt dt ^dr dt at dr' dr dr ^dcj) d<i) sin^ (j)de de'' 

where A", A*"*, A'"'" and A^ are smooth functions on the 2-jets of the bundle {r,A,B,C,D) -^ 
(r). In other words, of the ten components in the field equations, the dynamic reduction diagram 
automatically implies that 6 of these components vanish. Moreover, the reduced operator A is 
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constrained by the reduced Bianchi identities. Since dt and dr provide a basis for the invariant one 
forms on M, we know that the reduction of Divg S" is a hnear combination of dt and dr, 

Divg S = St dt + Sr dr. 

By direct computation, one finds that the dt and dr components of the reduced Bianchi identities 
are 



and 



i-^[7(2ylA'-* + i?A-)]=0 



-(-^[7(2CA'''' + BA''*)1 - iA" - C'A'^'' - BA''* - 2I)A") = 
2 ^ 7 dr ' 



where 7 = D\ —B"^ — AC. It foUows from the first of these identities and the transformation 
properties of A, B, A*"* and A*"^ under the residual scaUng t ^ et that 

2AA''* + BA'''' = 0. 

This same identity can be derived by first observing that the principle of symmetric criticality holds 
for the action G and then by applying Noether's second theorem to the reduced Lagrangian with 
symmetry ©off, 

Consequently of the four ODE arising in the stationary, spherically symmetric reduction of the 
field equations one need only solve the two equations 

A" = and A'''' = 0. 

The remaining two equations 

A''* = and A^^ = 

will automatically be satisfied (assuming D ^ 0, A ^ 0). We stress that these conclusions actually 
hold true for the stationary, rotationally invariant reductions of any generally covariant metric field 
equations derivable from a variational principle. | 

Example 6.6. Static, Spherically Symmetric Reductions. A metric is static and spherically 
symmetric if, in addition to being invariant under time translations and rotations, it is invariant 
under time reflection. The symmetry group G now includes the transformations t ^ t + c and 
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t ^ —t and therefore the isotropy subgroup Gx„ of the point xq = (to, xo) now includes the reflection 
t -^ 2to — t. The fibers of the kinematic bundle are now constrained by (6.22) along with 

bgb^ — 9, where fe = diag[— 1, 1, 1, 1]. 

This forces i? = in (6.23) so that the fibers of the kinematic bundle are now 3 dimensional and 
the general invariant section is 

ds^ = A{r)dt^ + C{r)dr'^ + D{r)dn'^. 

The automorphism group for this bundle is now r -^ f{r) and t ^ et and the A{M) summand 
in (6.27) does not appear. This example shows that while discrete symmetries will never result in 
a reduction of the dimension of the orbit space M, that is, the number of independent variables, 
discrete symmetries can reduce the fiber dimension of the kinematic bundle, that is, the number of 
dependent variables. 

Example 6.7. Plane Waves. As our next example from general relativity, we consider a class 
of plane wave metrics [12]. We take M ~ R^ with coordinates {u,v,x,y) and let P{u) and Q{u) be 
arbitrary smooth functions satisfying P'{u) > and Q'{u) > 0. The symmetry group on M is the 
five-parameter transformation group 

u' ^ u, v' = V + ei+SiX + e5y + l/2[s2S4: + £3S5 + £lP{u) + slQ{u)), (6.29) 

x' = X + S2 + £4:Piu), y' = y + ez + e^Q{u), 

d d d 

with infinitesimal generators Vi — — — , V2 — -r— , V3, — — — , 

ov ox oy 

V4^x—+P{u)— and V5^yTr + Qiu)Tr- 
ov ox ov oy 

The only non- vanishing brackets are 

[V2,V4]^Vi and [V3,V5] = Vi 

so that, regardless of the choice of functions P and Q, the abstract Lie algebras or groups are the 
same although the actions are generically different for different choices of P and Q. The coordinate 
function u is the only invariant and the orbits of this action are 3-dimensional. Therefore, at each 
point the isotropy subgroup is two dimensional and it is easily seen that, at xq = {uq, vo,xo,yo), the 
infinitesimal isotropy Fxo is generated by 

Zi = V4-xoVi~P{uo)V2 and Z2 = V5 - yoVi ^ Q{uo)V3. 
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At X G M the metric components g — [gij] of a G invariant metric satisfy the isotropy conditions 



gai + a\g — and ga2 + a\g — 0, 



where 



oi = 



0' 

10 

P'{u) 

0. 



and 02 = 



We find that the solutions to (6.30) are 
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P'{u) 







nu) 



Thus the kinematic reduction diagram is 



{u,A,B) 



{u) 



({M 



{x\A,B) 



ix^) 



id 



1 



and the inchision map l sends (A, B) to ds^ =- Adu^ + Bj, where 

„ , , dx'^ dy'^ 

j = -2dudv + —— + — 



P'{u) Q'{u) 
The most general G invariant metric is 

ds^ = A{u)du'^ + B(u)j. 



(6.30) 



(6.31) 



From the form of the most general G invariant symmetric type („) tensor, we are assured that the 
reduced field equations take the form 

ov ov ou ov ov ou ox ox ay ay 

Every G invariant one-form is a multiple of du so that there is only one non-trivial component to 
the contracted Bianchi identities and, indeed, by direct computation we find that 



DivsA^ = — (BA'^)dM. 
du^ ' 
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Since this must vanish idcnticaUy, wc conclude that the A'' component of the reduced field equations 
is of the form 

where c is a constant. Either the constant c is non-zero, in which case the reduced equations are 
inconsistent and there are no G invariant solutions, or else c ~ and the reduced equations consist 
of just the single equation A^^ — 0. For generally covariant metric theories the case c 7^ can only 
arise when the field equations contain a cosmological term [37]. 

It is easy to check that while the isotropy algebras Fxg are all two-dimensional abelian subalgebras, 
on disjoint orbits none are conjugate under the adjoint action of G. Hence the group action (6.29) is 
not simple and consequently the kinematic bundle for this action need not decompose according to 
(6.27). Indeed, the tensor 7 cannot be identified with any G invariant quadratic form on the orbits 

G/Gxq- I 

Example 6.8. Symplectic Reduction and Group Invariant Solutions. It is important to 
recognize the fundamental differences between symplectic reduction and Lie symmetry reduction for 
group invariant solutions of a Hamiltonian system with symmetry. Let M be an even dimensional 
manifold with symplectic form lo and let _ff : A^f ^ R be the Hamiltonian for a dynamical system on 
M . For the purposes of this example, it suffices to consider reduction by a one dimensional group 
of Hamiltonian symmetries generated by a vector field V with associated momentum map J, 

V-^Lo^dJ. (6.32) 

In symplectic reduction the reduced space M is obtained by [i] restricting to the submanifold of 
M defined by 

J = ^i = constant, 

and then [ii] quotienting this submanifold by the action of the transformation group generated by 
V. Both Lu and H descend to M and the reduced equations are the associated Hamiltonian system 
on M. Since diniM ~ diniAf — 2, the reduction in the number of dependent variables is 2. The 
solution to the original Hamiltonian equations are obtained from that of the reduced Hamiltonian 
equations by quadratures. 

To compare with symmetry reduction for group invariant solutions, we transcribe Hamilton's 
equations into the operator-theoretic setting used to construct the kinematic and dynamic reduc- 
tion diagrams. Let E — M x H ^ R he extended phase space so that the differential operator 
characterizing the canonical equations is the one- form valued operator on J^ (E) defined by 

A = X-'oj-dH. 
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Here X is the total derivative operator given, in standard canonical coordinates {u^,pi) on M, by 

_ d _ d .j9 . d 

dt dt du^ dpi 

It is not difficult to show that if V is any vector field on M, then the prolongation of V to J^{E) 
satisfies [X, pr^ V] — and therefore V^ is a symmetry of the operator A whenever V^ is a symmetry 
of Lo and H. 

Since 1/ is a vertical vector field on E it is "all isotropy" and the kinematic bundle is the fixed 
point set for the flow of V, 

The dimension of Hr{E) therefore depends upon the choice of V and is generally less than the 
dimension of E by more than 2 (the decrease in the dimension in the case of symplectic reduction). 
In short, it is not possible to identify the fibers of the kinematic bundle with the reduced phase space 
M . Moreover, from (6.32), it follows that points in Kr{E) always correspond to points on the singular 
level sets of the momentum map and, typically, to points where the level sets fail to be a manifold. 
Thus the invariant solutions are problematic from the viewpoint of symplectic reduction and are 
subject to special treatment. See, for example, [5] and [20]. Finally, there is no guarantee that 
the reduced equations for the group invariant solutions possess any natural inherited Hamiltonian 
formulation. 

We illustrate these general observations with some specific examples. First, if T^ is a translation 
symmetry of a mechanical system, then J is a linear function and symplectic reduction yields all the 
solutions to Hamilton's equations with a given fixed value for the first integral J. Since the vector 
field V never vanishes, the kinematic bundle is empty and there are no group invariant solutions. 

Second, for the classical 3-dimensional central force problem 

" = -I{p)u, V = -f{p)v, w = -f{p)w, 



where p = \/v? + w^ + w^, the extended phase space i<^ is R x R^ -^ R with coordinates 

{t,U,V,W,Pu,Pv,Pw) -^ (t), 

the symplectic structure on phase space is lo — du A dpu + dv A dpy + dw A dpw and the Hamiltonian 
is iJ = ^(pI +pI +pI,) + (j>(p), where (j)'{p) = pf{p). The vector field 

^^ d d d d 

^ "" ^"^" + ""5 P"^ ^^^"S — 

ov ou opv apu 
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is a Hamiltonian symmetry. 

The kinematic bundle for the V invariant sections of E is 

{t,W,Pw) ^ {t,W,Pw) > {t,U,V,W,Pu,Pv,Pw) 

TT TT TT (D.OOJ 

it) .-ii^ (t) -i^ (i), 

where i(i, w,pw) — (t, 0, 0, w, 0, 0, p^), the invariant sections are of the form 

t^(0,0,u;(t),0,0,p^(t)), 

and the reduced differential operator for the V invariant solutions is 

A = (w -pw)dpw - [pw +wf{\w\))dw. 

Let us compare this state of affairs with that obtained by symplectic reduction based upon 
the Hamiltonian vector field V. The momentum map associated to this symmetry is the angular 
momentum 

J = -upv + vpu. 

The level sets J = fi are manifolds except for fi = 0. The level set J = is the product of a plane 
and a cone whose vertex is precisely the fiber of the kinematic bundle. To implement the symplectic 
reduction, we introduce canonical cylindrical coordinates {r,9,w,Pr,P9,Pw), where 

u^rcosO, v~rs[n9, 

a PO ■ a ■ Q , P^ a 
Pu = Pr cos a sm tl , Pv = Pr sm -\ cos . 

r r 

Note that this change of coordinates fails precisely at points of the kinematic bundle. In terms of 

these phase space coordinates, the symplectic structure is still in canonical form lo = dr A dpr + d6 A 

dpe + dw A dpw, the Hamiltonian is H = —{p^ H — ttPq + p^) + (\>{\J r"^ + w^), and the momentum 

map is J = —pQ. We can therefore describe the symplectic reduction of E by the diagram 

{t,r,W,Pr,Pw) < {t,r,9,W,Pr,Pw) > {t,r,9,W,Pr,P0,Pw) 



it) < — (t) — > (t). 

The reduced symplectic structure is then uj ^ dr A dpr Adr + dw A dp^ , the reduced Hamiltonian is 
^ I „2 

H ~ —{Pr H — T +P11.) + Myr'^ + w^), and the reduced equations of motion are 
2 r^ 

„2 

r^Pr, Pr ^ -rf{\/r^ +w^) + —, w^pw, Pw ^ -wf{\/r^ + w"^). 

Given a choice of /i and solutions to these reduced equations, we get a solution to the full equations 
via 9 = —/it + const. | 
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7. Appendix. Wc summarize a few teclinical points concerning group actions on fiber bundles and 
the construction of the kinematic and dynamic reduction diagrams. For details, see [3]. 

A. Transversality and Regularity. Let G be a finite dimensional Lie group acting projectably 
on a bundle tt: E ^ M. We say that G acts transversally on E if, for each fixed p G E and each 
fixed g (z G, the equation 

""(ff ' p) — ""(j*) implies that g ■ p — p. (7-1) 

Thus each orbit of G intersects each fiber of E exactly once. For transverse group actions the 
orbits of G in i<^ are diffeomorphic to the orbits of G in M under the projection map t:: E ^ M. 
Projectable, transverse actions always satisfy the infinitesimal transversality condition (2.4) but the 
converse is easily seen to be false. 

Let us say that the action of G on M is regular if the quotient space M = M/G is a smooth 
manifold and the quotient map q^ '■ M ^ M defines M as a bundle over M. The construction of the 
orbit manifold M is discussed in various texts, for example, [1], [4], [29], [31]. The assumption that 
the action of G on M is regular is a standard assumption in Lie symmetry reduction. For simplicity 
we suppose that Af is a manifold without boundary but, as Example 6.4 shows, this assumption can 
be relaxed in applications. 

The fundamental properties of transverse group actions are described in the following theorem 
which is proved in [3]. 

Theorem 7.1 (The Regularity Theorem for Transverse Group Actions). Let G be a 
Lie group which acts projectably and transversally on the bundle n: E ^ M. Suppose that G acts 
regularly on M . 

[i] Then G acts regularly on E and E = E/G is a bundle over M . 

[ii] // the orbit manifold M is Hausdorff, then the orbit manifold E is also Hausdorff. 

[iii] The bundle E can be identified with the pullback of the bundle n: E ^ M via the quotient map 
(\m--M ^M. 

[iv] Let U be an open set in M and let U = q~/(t/). There is a one-to-one correspondence between 
the smooth G invariant sections of E over U and the sections of E over U. 

B. Transversality and the Kinematic Bundle. Lemma 3.2 implies that the action of G on 
E always restricts to a transverse action on the set Kg{E). In fact, it is not difficult to characterize 
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Kg{E) as the largest subset of E on which G acts transversally or, alternatively, as the smallest 
set through which all locally defined invariant sections factor. For Lie symmetry reduction without 
transversality the assumption that Kq (E) is an imbedded subbundle of E now replaces the infinitesi- 
mal transversality condition (2.4) as the underlying hypothesis for the action of G on E (together, 
of course, with the regularity of the action of G on M). In particular, the assumption that the 
dimension of Kq,^{E) is constant as x varies over M is clearly a necessary condition if one hopes to 
parameterize the space of G invariant local sections of E in terms of a fixed number of arbitrary 
functions. There are a variety of general results which one can apply to check whether Kg(E) is 
subbundle of E. To begin with, if x,y G M lie on the same G orbit, that is, ii y ~ g ■ x for some 
g (z G, then it is not difficult to prove that 

By virtue of this observation it suffices to check that the restrictions of Ka (E) to the cross-sections 
of the action of G on M are subbundles. For Lie group actions G which admit slices on M, it is not 
difficult to establish (see [4]) that the kinematic bundles for the induced actions on tensor bundles of 
M always exist. For compact groups acting by isometrics on hermitian vector bundles the existence 
of the kinematic bundle is established in [11]. 

Granted that Ka{E) — > M is a bundle. Theorem 3.3 now follows from Theorem 7.1. Theorem 7.1 
also shows that there is considerable redundancy in the hypothesis of Theorem 3.3. 

We emphasize that the action of G on i<^ itself need not be regular in order to construct a smooth 
kinematic reduction diagram. This is well illustrated by Example 19 in Lawson [26] (p. 23). 

C. The Bundle of Invariant Jets. The following theorem summarizes the key properties of 
the bundle Inv''(^) ^ M. 

Theorem 7.2. Let G be a projectable group action on t:: E ^ M and suppose that E admits a 
smooth kinematic reduction diagram (3.5). 

[i] Then Inv [E) is a G invariant embedded submanifold of J^{E). 

[ii] The action of G on Inv [E) is transverse and regular. 

[iii] The quotient manifold Inv {E)/G is diffeomorphic to j''{ka{E)) and the diagram 

J^{Rg{E)) ^^^^^ Inv'=(S) — '-^ J^{E) 

[ 1 1 
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commutes. 

This theorem impUes that the same hypothesis on the action of G on the bundle t:: E ^ M needed 
to insure that the kinematic reduction diagram is a diagram of smooth manifolds and maps also 
insures that the bottom row of the dynamical reduction diagram (4.2) exists. Therefore to guarantee 
the smoothness of the entire dynamic reduction diagram one need only assume, in addition, that 
I?i„v is a subbundle of T>. 

D. The Automorphism Group of the Kinematic Bundle. For computations of the auto- 
morphism group of the kinematic bundle it is often advantageous to use the fact that ©* fixes every 
G invariant section of E, that preserves the space of G invariant sections and that, conversely, 
under very mild assumptions, these properties characterize these groups. 

Theorem 7.3. Assume that there is a G invariant section through each point of Ko{E). Then the 
group 25* coincides with the subgroup of which fixes every invariant section of E, 

<&*={a£^\a-s — s for all G invariant sections s: M ^ E} 

and the group 25 coincides with the subgroup of which preserves the set of G invariant sections of 
E, 

= {aG©| a ■ s is G invariant for all G invariant sections s: M ^ E}. 
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